Hole pairing from attraction of opposite chirality spin vortices: 
Non-BCS superconductivity in Underdoped Cuprates 
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Within a gauge approach to the t — J model, we propose a new, non-BCS mechanism of supercon- 
ductivity (SC) for underdoped cuprates. We implement the no-double occupancy constraint with a 
(semionic) slave-particle formalism. The dopant in the t — J model description generates a vortex- 
like quantum distortion of the antiferromagnetic (AF) background centered on the empty sites, 
with opposite chirality for cores on the two Neel sublattices. Empty sites are described in terms of 
spinless fermionic holons and the long-range attraction between spin vortices on two opposite Neel 
sublattices serves as the holon pairing force, leading eventually to SC. The spin fluctuations are 
described by bosonic spinous with a gap generated by scattering on spin vortices. Due to the no- 
double occupation constraint, there is a gauge attraction between holon and spinon, binding them 
into a physical hole. Through gauge interaction the spin vortex attraction induces the formation 
of spin-singlet (RVB) spin pairs by lowering the spinon gap, due to the appearance of spin-vortex 
dipoles. Lowering the temperature, the proposed approach anticipates two crossover temperatures 
as precursors of the SC transition: at the higher crossover a finite density of incoherent holon pairs 
are formed, leading to reduction of the hole spectral weight, while at the lower crossover a finite 
density of incoherent spinon RVB pairs are also formed, giving rise to a gas of incoherent preformed 
hole pairs with magnetic vortices appearing in the plasma phase, supporting a Nernst signal. Finally, 
at an even lower temperature the hole pairs become coherent, the magnetic vortices become dilute 
and SC appears beyond a critical doping. The proposed SC mechanism is not of the BCS-type, 
because it involves a gain in kinetic energy, due to the lowering of spinon gap, and it is "almost" 
of the classical 3D XY-type. Since both the spinon gap, describing short-range AF order, and the 
holon pairing, generating SC, originate from the same term in the slave-particle representation of 
the t — J model, the proposed approach incorporates a strong interplay between AF and SC, giving 
rise to a universal relation between Tc and the energy of the resonance mode, as observed in neutron 
scattering experiments. 

PACS numbers: 71.10.Hf, 11.15.-q, 71.27.-|-a 



I. INTRODUCTION 

The high temperature superconductivity (SC) in 
cuprates, discovered 25 years agoji still remains a major 
challenge in the condensed matter physics. In spite of the 
enormous progress made in materials synthesis, crystal 
growth, experimental studies of physical properties and 
theoretical interpretation, there is still no consensus yet 
regarding the anomalous normal-state properties and SC 
mechanisms in these cuprate compounds. There is a re- 
cent review article^ on various approaches attacking this 
extremely difficult problem, including Resonance Valence 
Bond (RVB) slave particle gauge approaches, spin fluc- 
tuation models, stripe models, phonons, three-band sce- 
nario, etc. We share many viewpoints expressed there. 



and to save space, we refer the readers to that review 
article, not repeating those comments here. 

In this paper we propose a new mechanism of SC in 
hole-underdoped High Tc cuprates, using the spin-charge 
gauge approach to the 2D t-J model, describing the Cu- 
O planes.'^ In this approach the t-J model (with t/J 
as the main parameter) satisfying the single-occupancy 
constraint, is treated systematically within the same set 
of approximations, to study both normal-state and SC 
properties. The exchange J-term giving rise to antifer- 
romagnetism (AF) is also serving as the "glue" leading 
to SC, thus implementing the interplay of AF and SC 
in an explicit form. The proposed SC mechanism is not 
of the BCS-type, and it involves a gain in kinetic en- 
ergy by lowering the spinon gap due to the appearance of 
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spin-vortcx dipolcs. The main features of this non-BCS 
description of SC are consistent with the experimental 
results in underdoped cuprates, including a natural d- 
wave SC pairing parameter, and especially the contour 
plot of the Nernst signal<^i^ We can also derive the SC 
transition as "almost" of the classical 3D XY-type, while 
the calculated transition temperature shows a universal 
ratio to the resonance mode energy observed in neutron 
experimentsi^ 

Our formalism basically belongs to the "strong 
correlation-slave particle tent", where, a U(l) field is in- 
troduced to gauge the global charge, while a SU(2) field is 
introduced to gauge the global spin. Through the gauge 
field, a vortcx-likc quantum distortion of the AF back- 
ground is generated around the empty site (described in 
terms of fermionic spinless holon) with opposite chiral- 
ity for cores on two Neel sublattices. In the presence 
of these vortices the spin excitations (bosonic spin 1/2 
spinon), originally gapless without doping, correspond- 
ing to a long range AF order, acquire a finite gap due 
to scattering on these vortices (similar to the wave lo- 
calization of light propagating in random media), con- 
verting LR AFO to a short range order (SRO). Within 
this approach the physical hole is a bound state of holon 
and spinon with a "glue" (binding force) coming from 
an emergent U(l) slave particle gauge field. Here the 
spinon and holon are neither confined (as in the ordi- 
nary Fermi liquid), nor decoupled (as in ID t~J model 
in the small J limit), but rather forming a "composite 
particle" -physical hole in a strongly correlated system. 
It is not anymore a "neat" quasiparticlc, but rather hav- 
ing a strongly temperature-dependent lifetime due to the 
gauge field (coupled to holons with a finite Fermi sur- 
face). Similarly, the magnon is a composite particle made 
of spinon and antispinon, again with the "gauge glue" . 
In fact, the "composite" characteristics are responsible 
for all exotic properties in the "pseudogap phase" (PG), 
in particular, the interplay of the SR AFO (exhibited 
as a finite magnon mass gap) with the dissipativc mo- 
tion of charge carriers, showing up as lifetime effect of 
the physical hole, results in a metal-insulator crossover, 
a pronounced phenomenon in the underdoped cuprates. 
A number of peculiar features of cuprates in the normal 
state can be well explained within this schemcj^ Here this 
approach is generalized to consider the SC state. 

The gluing force of the SC mechanism is an attraction 
between holons generated by spin vortices on two op- 
posite Neel sublattices, centered around the empty sites 
(holes). This attraction which shares the same origin of 
spin exchange J-tcrm leading to AFO, was neglected as a 
subleading term in considering the normal state proper- 
ties. Physically, the hole is assigned an additional "pscu- 
dospin" index marking the belonging Nccl sublatticc in 
a range characterized typically by the AF correlation 
length, as long as the SR AFO persists. This attraction 
describes the tendency towards vortex-antivortex bind- 
ing, or reduction of the AF exchange energy loss. In fact, 
the formation of the vortex dipolcs effectively reduces the 



density of free vortices scattering the spin waves, and as 
a consequence the kinetic energy of spinon gains and the 
AF correlation is enhanced as well. To materialize the SC 
transition we propose the following three-step scenario: 

At the highest crossover temperature, denoted as Tph, 
a finite density of incoherent holon pairs are formed. We 
propose to identify that temperature with the experi- 
mentally observed (upper) PG temperature, where the 
in-plane resistivity deviates from the linear behavior. A 
BCS-like d-wave pairing of holons is derived by "super- 
posing" two p-wave like pairing in a reduced Brillouin 
zone. However, the holon pairing alone is not enough for 
SC to appear. Again, through the "gauge glue" coming 
from the U(l) slave particle gauge field, so crucial for 
the interpretation of exotic properties of the PC phase, 
the spin vortex attraction induces the formation of spin- 
singlet (RVB) spinon pairs with a reduction of the spinon 
gap. Physically, the spinous will feel "less disturbance" 
from the formation of vortex-antivortex pairs (dipoles) . 

At the intermediate crossover temperature, denoted as 
Tps, a finite density of incoherent spinon RVB pairs are 
formed, which, combined with the holon pairs, gives rise 
to a gas of incoherent preformed hole pairs. We pro- 
pose to identify that temperature with the experimental 
crossover corresponding to the appearance of the Nernst 
signal. The calculated contour plot of the spinon pairing 
parameter is compared with the Nernst signal contour 
plot,— showing a good agreement. 

Finally, at an even lower temperature, the SC transi- 
tion temperature Tc, both holon pairs and RVB pairs, 
hence also the hole pairs, become coherent. It will be 
shown that the phase coherence is established via a phase 
transition of a planar gauged quantum XY-type, almost 
identical to that of the classical 3D XY-model. The 
SC transition temperature is calculated as a function of 
doping concentration, and is compared with the scaled 
value of the resonance mode energy, observed in neutron 
experiments^ to show the universal ratio between these 
two quantities anticipated from our theoretical treat- 
ment. 

The rest of the paper is organized as follows. Section II 
is a brief introduction to our scmionic spin-charge gauge 
approach, to make the paper more self-contained. Sec- 
tion HI is devoted to the holon pairing mechanism. In 
section IV we discuss the spinon pairing, while in section 
V the SC transition is considered. Discussions and con- 
clusions are given in section VI. Several technical deriva- 
tions are outlined in Appendices. A preliminary report 
of the present work has already appeared in Ref. 0- 



II. THE SPIN-CHARGE GAUGE APPROACH 

A. Slave semions 

In this subsection we present an outline of our slave- 
particle approach using the "semionic" spin-charge de- 
composition, applicable only in 2D (and ID) systems. We 



3 



assume that the main features of the low-energy physics 
of the hole-doped cuprates can be captured by the t — J 
model 

+ JY,iS^■S,+lhS,)]PG, (1) 

where Pq is the Gutzwiller projection imposing no- 
double-occupation condition and the lattice sites corre- 
spond to those of the Cu atoms in the Cu02 planes of 
the cuprates. The particle number and spin operators 
are defined as 

^i=Yl ^t^^'^^ = ^L-^c,;3 ■ (2) 

t, J, and /i in Eq. ([l} are hopping amplitude, spin ex- 
change and chemical potential, respectively. The hole 
operator carries both charge and spin degrees of free- 
dom, with no-double-occupation constraint. Formally 
they can be treated separately by the standard slave- 
particle approach, Ci„ = hi„h\^ where hi is a fermionic 
holon operator, and bia is a bosonic spinon operator. 
The no-double occupation condition is automatically en- 
sured by the spinless fermion, while the correct counting 
of degrees of freedom (dof) is imposed by the constraint 
^L^io- = 1 on the b field, so that Ui = hih\ = l — h\hi. 
At half filling the charge degree of freedom is frozen, then 
the slave-particle transformation reduces to the standard 
Schwinger-boson approach. 

In 2-1-1 dimensional systems, one can bind statistical 
fluxes to particle-excitations, resulting only in a change 
of the statistics. This is achieved in the Hamiltonian for- 
malism by minimally coupling the matter fields to suit- 
able composite "statistical gauge field operators". The 
introduction of these fluxes in the lagrangian formalism 
is materialized via statistical Chern-Simons gauge fields. 
In our case, holes carry both charge and spin degrees of 
freedoms, so we associate two statistical gauge fields with 
hole operators, one of which is a U(l) gauge field B cou- 
pled to the holon /i-field and related to the charge, while 
the other is an SU(2) field V coupled to the spinon b- 
ficlds and related to the spin. By carefully choosing the 
coupling constants of the corresponding Chern-Simons 
terms, we can keep the original hole field still fermionic. 
In the Hamiltonian formalism, the statistical gauge fields 
operators can be chosen as follows: 

Bt,{x) = i^fi;-€'^arg(f - 0, 

r 

^ i 

where the sums are carried over lattice sites I, the sum 
over the spin indices (a, b = x^y, z) is understood here- 
inafter and the function arg(af) is the angle of the vector 



X. The corresponding U(l) and SU(2) fluxes, and ^s, 
bound to the hole at site j are given by 

{e''''^'^U^{Pe'^r''-'^^(^\p 

= ie^^rS^r^'-^(^-^U, (4) 

where a,/3 = 1,2 are the SU(2) spin indices. The in- 
tegration runs over a path joining j to infinity and P 
denotes the path-ordering. Binding the holon to the 
U(l)-fiux generated by B and the spinon to the SU(2)- 
fiux generated by V chosen as in Eq. ^ one obtains 
U(l) and SU(2) invariant fields, respectively, both obey- 
ing semionic statistics,^^ i.e. their interchange produces 
a ±i factor, intermediate between the bosonic -|-1 and 
the fermionic —1 case, whence the name "scmion"jiS 

This "semionic" approach is quite suitable to study the 
physics of holes dressed by a spin vortex as described in 
the Introduction because the SU(2)-gauge field naturally 
incorporates the spin vortices. To show that is indeed 
the gauge field associated with spin vortices let us con- 
sider the simplest case of one hole located at / with spin 
5°. Then V^{x) simplifies to S'°i9^t arg(a? — l) and using 
£f,^d^d^ arg(x - y) = 6{x - y) we get 

e^''(a^C(2^)) = ^"'5(f-r). (5) 

Eq. ([5]) is a spin analogue of the charged vortex intro- 
duced by Laughlin in the fractional quantum Hall effect 
and, in fact, a semionic representation of the electron 
was advocated originally by him in the early days of high 
temperature SC boomj^i 

B. Improved Mean Field Approximation 

In this subsection we sketch the key approximations 
involved in our approach to the "normal" state; one of 
these approximations (the optimization procedure) ap- 
pears rather unconventional in slave-particle approaches. 

Being too difficult to be solved exactly, the gauge field 
approach outlined above provides a reasonable base of 
an improved mean field analysis^ that, dimensionally re- 
duced, works quite well for one dimensional t—J modelfi^ 
reproducing correctly also the non-trivial critical expo- 
nents of its correlation functions (the spin-vortices be- 
come kink strings in ID). In two dimensions, this mean 
field theory involves an optimization of the spin config- 
uration around holons dressed by vortices, although it 
can be carried out only approximately and not rigor- 
ously as in the one- dimensional case. In the improved 
semionic mean field approximation (MFA), the spinon 
configurations around holons arc optimized leading to a 
new bosonic spinon on the optimized spinon-background, 
denoted by z, which is therefore different from the 6-field 
in the standard slave fermion approach, but still satis- 
fying the constraint zj^Zia = 1- From now on it is this 
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spinon that we refer to. In the adopted MFA we neglect 
the holon fluctuations in B and the spinon fluctuations 
in V. This leads to a much simpler form of the two sta- 
tistical gauge fields denoted by B and V, respectively, to 
distinguish them from the exact values. The B field is 
actually a static one, without dynamics, 

^t^^l^^t^^^Six-l), (6) 
/ 

and it provides a tt-Aux phase factor e*^'^ per plaquette 
for the holon field because Jq Bij = tt. For the SU(2) 
gauge field only the a^-component survives: 

V^i^) « - E f^Vhi^-^d^ arg(f - f), (7) 
I 

with a pure-gauge static term being gauged away. Note 
that there is the holon number operator in the right hand 
side of Eq. ([7]), which means that the spin vortex is always 
centered on the hole, and its topological charge (named 
chirality) is (—1)''' depending on the parity of the site 
index, where \l\ = + ly The effect of the optimal spin 
flux is then to attach a spin-vortex to the holon, with 
opposite chirality on the two Neel sublattices, and the 
rigidity holding up a vortex being provided by the AF 
background. These vortices take into account the long- 
range quantum distortion of the AF background caused 

I 



The Euclidean Lagrangian used in the path-integral 
formalism is then obtained by replacing the field oper- 
ators h,h^ and i,zt with Grassmann {h,h*) and com- 
plex number (z,z*), respectively, and adding the time- 
derivative terms 

Y^'hdoh, + h*h,){~lf\z*doz,. (11) 

i 

The Hamiltonian Eq. ^TU\i is our starting point for de- 
scribing the High Tc cuprate SC. At the mean field level, 
the first two terms describe the motion of the holons, 
which are coupled to the spinous through the AM fac- 
tor whose modulus we treat as a constant, giving a small 
correction to the hopping amplitude t of holons that we 
neglect. Its phase factor ^ e*'''^ instead cannot be ne- 
glected, and it provides a gluing force between the spinon 
and holon. Then the mean field Hamiltonian of holon 
reads 

H^. = tJ2hy^^'^+0'^^h.. + h.c^^^Y,h^h,. (12) 

(id) » 



by the insertion of a dopant hole, as first discussed in 
Ref. [13I As in the one-dimensional case the optimization 
involves also a spin-flip associated to every holon jump 
between different Neel sublattices, hence in the i-term the 
spinous appear in the "ferromagnetic" AfHeck-Marston 
(AM)ii form xtj = {zje'^^'^ Zj)*^'\ where #(i) denotes 
complex conjugation if i belongs to the "odd" sublattice, 
with a phase ambiguity left by the optimization, whereas 
in the J-term it appears in the "AF" RVB form Af^ = 

e"''z,„(e'^-i^'"-%)/3, where 

V!^ ^ [ dx^V^{x)^V^{^). (8) 

The above AM/RVB dichotomy is peculiar to the semion 
approach involving the SU(2) spin rotation group even in 
ID, where it can be rigorously derived. It does not appear 
in the standard U(l) slave fermion or boson approaches. 

In the above MFA the hole field operator can be decom- 
posed as a product of the holon and the spinon operators 
along with fluxes: 

c,. = /ile'*'(e'*'z.)^- (9) 

The resulting MFA of the J model Eq. ^ is written 
in terms of holon fields hi and spinon field Zi as 



(10) 

I 

In two-dimensional bipartite lattices for fermions in 
magnetic field the optimal flux per plaquette is tt at half- 
filling (Lieb theoremi^) and numerically it is also true for 
close fillings at low temperatures, whereas it is zero suffi- 
ciently far away from half-filling. Therefore the optimal 
fiux in a plaquette for {Bij -\- Oij) is arguably tt for small 
doping and low temperatures, and for sufficiently high 
dopings and/or high temperatures. We conjectured that 
this corresponds to the crossover between the "pseudo- 
gap phase" (PG) and the "strange metal phase" (SM) as 
varying the doping or temperature in the cuprates, where 
PG is the "lower pseudogap" in the literature identified 
with the inflection point in resistivity and the broad peak 
in the specific heat coefficient 7. This conjecture is sup- 
ported by the comparison of the behavior of the theo- 
retically derived crossover temperature T* « -^|ln(5|p 
with experimental data, where the appearance of |ln5| 
is due to the long-range tail of spin vortex interactions. 
Therefore we fix the phase ambiguity left by the opti- 
mization in the AM term by choosing this phase zero for 
PG since B has already tt flux and is opposite to B for 
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SM to effectively cancel B. 

If we replace the holon density by its average in MFA, 
the third term in Eq. (fTO|) describes the motion of z- 
spinons with J renormalized to J = J(l — 26). Without 
doping, using the identity 

|A^P + |X|,P = 1, (13) 

holding for bosonic spinous, together with Eq. in the 
continuum limit it leads to a standard nonlinear cr-model 
describing the low energy physics of the AF background. 
With doping the spinous are scattered by holons dressed 
by spin vortices and that leads to a short range AF corre- 
lation. Such a process is revealed by expanding the SU(2) 
phase factor inside the RVB factor in the third term of 
Eq. (jlOp to the second order, obtaining in the continuum 
limit, self-consistently in the region with unbroken SU(2) 
spin symmetry 

J J <fxV^\x)zl{x)z^{x). (14) 

In MFA we replace V^^, positive definite by definition, 
by a statistical average. The spatial average of V^^(x) at 
fixed holon position Xi by using Eq. ([7|) reads 

Y^yif\+\=\A-\x,^x,), (15) 

where |i| = \xi\ and A is the two dimensional lattice 
Laplacian.— Eq. (fT5|) appears as the energy of a two di- 
mensional Coulomb gas with the lattice spacing as an 
ultraviolet cutoff, which can be evaluated at fixed den- 
sity (5 by a quenched approximation leading to a doping- 
dependent mass term for spinon, which in the low doping 
limit is given by 

ml{6)^{V^^)^\\5\n6l (16) 

consistent with AF correlation length {£,af ^ ("t-s)"^) 
at small 6 extracted from the neutron experiments In 
^AF the factor is just the mean distance between 
holes, while the factor |ln5| comes from the long-range 
tail of the vortex interactions and it turns out to be a 
key feature in many physical quantities in our approach. 
The spinon gap is also crucial for eliminating the over- 
counting of low-energy degrees of freedom often encoun- 
tered in slave-particle approaches, giving rise to prob- 
lems in the computation of thermodynamic quantitiesJ^ 
In fact, because of the spinon gap, the low-T thermody- 
namics in our approach is essentially dominated by the 
gapless holons, while the contributions of the transverse 
and scalar gauge fiuctuations to the free energy almost 
cancel each other In the Lagrangian form, our massive 
cr-model derived from Eqs. PT|) . dH]) and (|TE| can be 
conveniently written as 

J^s ^ - [ d^x[\{do-iAo)z^\'' -vl\{d^-iA^)z^\'' 



+ml{5)zlz^]{x), (17) 

where an implicit momentum cutoff is implied inside the 
magnetic Brillouin zone (MBZ), g = SJa^ with a lattice 
spacing, Vs — Ja, and the emergent gauge field is 
generated by the fluctuations of spinous: 

A^^e'Q-^^zi{x)d^z^{x) + ..., (18) 

with Q the AF wave vector, and it corresponds to the 
long wavelength limit of 9ij , the phase factor of AM fac- 
tor Xij- Note that in the massive cr-model Eq. (fT7| . the 
constraint z^ z = 1 on the z-field is relaxed. Holons and 
spinous are coupled by the gauge field A^, giving rise 
to overdamped resonances for holes and magnons with 
strongly T-dependent life-timei^ This dependence origi- 
nates from the dynamics of the transverse mode of the 
gauge field that is dominated by the contribution of the 
gapless holons. Their Fermi surface generates an anoma- 
lous skin effect, with momentum scale 

Q « (t4)1/3, (19) 

known as the Reizer momentum)^— where kp is the 
holon Fermi momentum measured from the Dirac point 
in TT-ffux phase. For the appearance of Reizer skin effect 
the presence of a gap for spinous is crucial, because gap- 
less spinous would Bose condense at low T thus gapping 
the gauge field through the Andcrson-Higgs mechanism 
and destroying the T-depcndent skin effect that reduces 
the coherence of hole and magnon. The transport physics 
of PG is dominated by the interplay between the short- 
range AF order due to spinous and the thermal diffu- 
sion induced by the gauge fluctuations triggered by the 
Reizer momentum, producing in turn the metal-insulator 
crossover!^ 

More generally, the above semionic mean field treat- 
ment based upon a spin-charge gauge approach to the 
t — J model provides a description of many transport 
and thermodynamic properties of High Tc cuprates in 
PG region^ii^i^i"— whose doping-temperature behavior 
is in qualitative, and sometimes even semi-quantitative, 
agreement with experimental data. In the following, we 
present details of the novel non-BCS mechanism for high 
Tc SC outlined in the Introduction. Here we just rewrite 
the SC order parameter in the approximation adopted 
above: 

A?, = e„/3Q„c,/3 ^ hlh] Af^e^^^^ , (20) 

which can be obtained by Hubbard-Stratonovich trans- 
formation in the path integral formalism. In the next 
section we discuss the holon pairing {hlh]), while m sec- 
tion |TV] the spinon pairing (A| ). 
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III. HOLON PAIRING 
A. Holon Hamiltonian with Attractive Interaction 

The Hamiltonian Eq. (fT2|) in PG describes the mo- 
tion of liolons which are subjected to a staggered tt-Aux 
field and the gauge field Oij, coupling them to spinous. 
To get the low energy physics of holon, we first neglect 
the gauge field 9ij generated by spinous, and it will be 
reinserted (in an approximate form) by Peierls substitu- 
tion. The remaining terms can be solved exactly. We find 
that the holon spectrum involves two Dirac cones due to 
the presence of the tt-Aux (Hofstadter mechanism).— The 
Fermi surface of holon is a small one with Fermi wave- 
vector kp ~ 7r(5.— Due to the staggered tt-Aux, we divide 
the square lattice into two sublattices, A(even sites) and 
B(odd sites). On each sublattice, the holon's annihila- 
tion operators are denoted by a and 6, respectively. The 
Hamiltonian Eq. (|12p of free holon can then be recast in 
a quadratic form 

^o' - E(%«lA- + ^■^■) - mE(4«a; + ^A-)> (21) 
fc fc 

where the momentum runs within the magnetic Brillouin 
zone (MBZ) and t^: = 2t{cosk^e''^/'^ + coskyer''^/'^). It 

is straightforward to obtain the spectrum e(k) = ±|ij;|, 
with the Fermi surface consisting of four half circles 
around (±7r/2, ±7r/2), as shown in Fig. [1^, where the 
blue lines are the boundary of MBZ. The Fermi energy is 
t5 approximately. There are two primitive reciprocal vec- 
tors, 7r± = (±7r, tt) by which we can translate the MBZ 
in the 3rd and 4th quadrants to get another equivalent 
rectangular one as shown in Fig. [T}d, which consists of 
two Dirac cones centered around Ql = (— 7r/2, 7r/2) (left) 
and Qr = (7r/2,7r/2) (right), respectively. In this trans- 
formation, wc note that a,- , - = ar and hr , ^ = —for, 

' k-\-7r± k fc+7r± fc' 

where a minus sign appears for 5-field defined on odd 
sublattice, but the form of Hamiltonian Eq. (PT|) is still 
invariant, because also changes sign after translation. 

Accordingly, all the holon operators can be labeled by 
the "flavor" index a = L,R distinguishing left and right 
Dirac zones, so does the Hamiltonian of free holons Hq = 

Ec.=L,K^o,a with 

Hn„ = y"{t r + h.c.) 

^>^<^ / a,fc fc ' 

fc 

a,k a,k a,k a,k'^ V ' 

where tj^j: fa 2t[-'kx + iky) and j: ~ 2t{kx + iky). In 

Eq. ()22|) , the momentum k only takes values in the range 
[— 7r/2,7r/2] x [— 7r/2, 7r/2], which is one quarter of the 
original BZ. 

Now we consider the holon-holon interactions. As 
shown in Sec. Ill B[ the last term in Eq. (fTO|) is repul- 
sive for holons which cannot be the pairing force between 
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FIG. 1: The Brillouin zone and Fermi surface of free holon 
with TT flux. The folded MBZ with blue lines as boundary in 
(a) is equivalent to the rectangular in (b). 

holons, and is usually negligible in the low doping limit. 
Meanwhile, the third term in Eq. (jlOp implies an effec- 
tive long-range interaction between the holons mediated 
by the spin vortices bound to holons, which turns out 
to be attractive between different Ncel sublattices. This 
leads to the instability of holons towards pair formation 
and is our key attractive force. Such an effect in the 
simplest form was first realized by S. TrugmariS^ in the 
early days of High Tc research, who pointed out that 
putting two holes next to each other on two Ncel sublat- 
tices would save energy J. We include this effect in MFA 
by introducing a term coming from the average of z^z in 
(|T4l) obtaining: 

JeffE(-l)''''"''^"'(*-i)^l^^^l'*J- (23) 
id 

In the static approximation for holons, Eq. ([23|) describes 
a 2D lattice Coulomb gas with charges ±1 depending on 
the Neel sublattice and coupling constant Jeff ~ J{z'^z), 
where the average {z'^ z) can be estimated from the free 
spinon spectrum (which will be given in the next section, 
see Eq. ((50)) by setting Aq = 0) with the following result 

= J{l-25){^K^+ml~m,), (24) 

where A is the momentum cutoff for spinon excitations. 
For 2D Coulomb gases with the above parameters a pair- 
ing appears for a temperature Tph ~ Jes/2-K (a more 
precise estimation is given later, in fact, Tp^ is the up- 
per PG crossover temperature determined by l^^lkp) of 
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Eq. (|34)) ') which turns out to be inside the SM "phase". 
Hence the whole PG "phase" hes below Tpf^. However, 
we will discuss only the SC arising from the PG phase, 
anticipating that extrapolation to SM phase will intro- 
duce only quantitative changes (actually the role of a 
next nearest neighbor hopping t'-term appears relevant 
in SM22,). In the large-scale limit the two dimensional 
Coulomb interaction gives rise to a screening effect, with 
a screening length is which in the Thomas- Fermi approx- 
imation is proportional to 1/ ^/S^ In view of the above 
considerations we approximate the large-scale effective 
potential in momentum space by 



Vesip) 



J, 



cff 



-2 • 



(25) 



The large-scale holon interaction then has the following 
simplified form 



E 

piqiP2q2 



92 j 



Pi "qi "q2 ^P2 



(26) 



Due to the long range tail of vortex-vortex interac- 
tion, the pairing strength for large momentum((f ~ (tt, 0)) 
transfer between different Dirac cones is much smaller 
than that for small momentum(g ~ 0) transfer. Hence, 
in the presence of interaction, the left and right flavors of 
holons can still be decoupled approximately. Considering 
the BCS approximation, where pairing occurs between 
holons in states with opposite momentum, one obtains 
the decoupled Hamiltonians Hj^ for each flavor a, 



We shall now focus only on the quasiparticles near the 
Fermi circles, which allows us to make the following gauge 
transformations for the holon operators with different fla- 
vors separately 



a.k a.k ^ ' a,k a,k 



,£/2 



(28) 



where the angles 6^ j are chosen to cancel the phase of 



i J so that the kinetic term reads 



~ ''Fkial^b^ j: + h.c.) 



(29) 



with vp = 2t. Eqs. ([?7)) and (P^]) are our basic equations 
to describe the pairing of holons. 



B. D-wave Pairing 

In this subsection, we show that the d-wave pairing 
symmetry is composed naturally of two p-wave pairing 
within the left and right Dirac cones, an idea first pro- 
posed by Sushkov et al^'^ in a different setting. The 



corresponding pairing parameter has a form respecti 
the Ciy rotation symmetry. 



A" - 



A''(fc)^ 
A'>(fc)^ 



if a = i?, 
■, if a = L, 



(30) 



where the momentum k is measured from Qrx^ respec- 
tively, and the magnitude of the order parameter is the 
same for both R and L flavors. Note that we are now 
working with the rectangular magnetic Brillouin zone 
(see Fig[TJ)) and the p-wave pairing takes place within 
the two circular Fermi surfaces. If transformed back to 
diamond magnetic Brillouin zone as in Fig. [T^, the order 
parameters in region HI and IV change their signs due 



to the fact that for 



-b-r, which leads to a perfect 



/c+7r± A;' 

d-wave pairing in the original Brillouin zone. 

Applying the standard BCS treatment we get the fol- 
lowing MF Hamiltonian 



fjh frh 



Z — a.k a,k a, — k ' 



(31) 



where the order parameter satisfies the gap equations. 



A'* - 



^\k = E ^e.(fc - <f)^^ tanh [^) . (32) 



It turns out that Eq. ([31]) has two decoupled branches of 
solutions (see Appendix]^ for details). One of them with 
higher energy without FS provides a matrix element sup- 
pressing the spectral weight of the original holon field h 
outside the MBZ as in PGi^ The other one is responsible 
for the low energy physics of holon pairing which we will 
focus on in the following and its spectrum has a simple 
BCS form 



(33) 



As common for non- weakly coupled attractive Fermi sys- 
tems, the MF temperature at which A'' becomes non- 
vanishing should be identified with the pairing tempera- 
ture Tph. 

For brevity, we consider the p-wavc order parameter in 
the right cone, which has the form A'i = AQ(fc)(cos6'^ — 

sin 6*^). The radial part AQ(fc) is decoupled from its angu- 
lar part approximately (see Appendix |^ , which is plot- 
ted in Fig. [2] for different values of the screening length 
is- We observe that holons near the Fermi surface take 
part in pairing which results in a peak of A''(fc) centered 
around k ^ kp. Actually, the number of holons partic- 
ipating in pairing is determined by the screening length 
£s- If we increase is, a higher percentage of holons can 
interact with each others at longer distance, the peak 
of A'*(fc) in Fig. [2] becomes higher and wider, that im- 
plies a bigger fraction of holons is involved in pairing. 
A more rigorous treatment would actually involve taking 
into account self-consistcntly the UV cutoff and chemical 
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potential change, as discussed e.g. in Refs. 130113 ll . but for 
simplicity we refrain to do that, assuming that our sys- 
tem is sufficiently BCS-like and our treatment catches 
already the key behavior, as Fig. [5] suggests. 
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FIG. 2: Plots of pairing gap of holon as functions of momen- 
tum in the upper panel, and as functions of temperature in 
the lower panel for different screening length £s. Aq is plotted 
in units of i. It is seen that the holons near the Fermi surface 
take part in pairing, leading to a peak of A(fc) center around 
k ~ kp. 

The maximum value of the order parameter at zero 
temperature can be taken as the typical energy scale of 
pairing strength, which is the value of Aq at the Fermi 
momentum. Though it is difficult to get an analytical 
solution of AQ{kp) from the radial gap equation (see 
Eq. (jA9p ). one can get an approximate expression for 
it as a function of the parameters Jeff, kp and i^, which 
has the following form 



A^(fcF)«0.06Jeff(fcF^s)exp 



(34) 



being not sensitive to the energy cutoff as long as the 
screening length £s is larger than 1/A. 

Now we can write down the d-wave order parameters 
near the original four Fermi arcs (see Fig. [T^) 

• « v\{kx — ky)/^/2 in quadrant I, 

• A^ « v^{—kx — ky)/^/2 in quadrant II, 



{—kx + ky)/^/2 in quadrant III, 



• A^ 



v\{kx -f ky)/\/2 in quadrant IV, 



where fA = \/2AQ{kp)/kp. 

So far we discussed the d-wave paring symmetry in the 
momentum space in the long wavelength limit, and now 
we check that when extrapolating the result to the lat- 
tice scale we recover the desired pairing symmetry in real 
space. Computing the nearest neighbor pairing between 
site X and x + 5 we get. 



1 



k-S 



(35) 



where V is the volume of the system and the summa- 
tion over k is in the range [— 7r/2,7r/2] x [— 7r/2, 7r/2]. 
Note that {ba.-q'aa,q) has the same symmetry as A^ - 
(sec Eq. (|A4[) '). then by using Eq. (pUj) . one can easily 
prove (bsag^g) = {bsa,-g_g) and {bsas+Si) = -{bgas+e^)^ 
which are the typical features of d-wave order parameters 
in real space. 



C. Nodal approximation and Gauge Invariance 

In the BCS approximation discussed in the previous 
subsection the holon is gapless only at the 4 nodal points 
of A|. However, in a large-scale gauge- invariant treat- 
ment whereas one can keep the modulus of the order 
parameter A'' as in BCS, we must include its spatially 
dependent phase, which we denote by (j>^{x). (A pre- 
cise procedure to go from the lattice to the continuum 
phase field is discussed in Ref. HI.) The effects of (j)'^{x) 
on holons is non-trivial and will be discussed in detail 
in Rcf. [261. However, to derive our basic RVB gap equa- 
tion in the next section it turns out that we can assume 
consistently that (f)^{x) doesn't break the nodal struc- 
ture. In fact the nodal structure appears if we neglect 
the phase fluctuations, so that the holon pairs are as- 
sumed condensed. According to Refs. [33l - l35l this is the 
correct procedure to deal with the gap equation for the 
modulus of the order parameter. However, if holon pairs 
are only formed but not yet condensed, it is incorrect to 
identify A'* as the gap for holons (see Ref. [26l ). In this 
subsection we utilize the Peierls substitution to take the 
gauge flelds back into account around the nodal points, in 
agreement with the above remarks. In the nodal approx- 
imation the momenta are expanded around the nodes in 
the four quadrants of the MBZ. In Fig. [31 we plot the 
nodal coordinate system in the 1st quadrant, where 



k. 



V2 



V2 



(36) 



In terms of fc+ and fc_ , using the the gap dependence on 
momentum k obtained in the last subsection, the energy 
spectrum around the node of the 1st quadrant is simply 
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FIG. 3: The coordinate system (k+,k-) taken in the nodal 
approximation. This is for the first quadrant. 



±yjvpk'^ + w^fc^, which arises from the nodal Hamilto- 
nian in the 1st quadrant 

^Ist nodal = VEk+f^z + ^'Afc-C^y, (37) 

which reproduces the spectrum of the gapless nodal ex- 
citation. Therefore, adding also the contribution of the 
phase (t)^{x) of the order parameter, a large-scale h/s 
gauge invariant Hamiltonian in real space reads 

n,h _(-id+-A+ + AQ -e'^'d- \ 

e-*^9_ ^^+-A^-A,)' ^^^^ 

where the emergent gauge field A^ and Ao(scc Eq. (fT8|) ) 
is reinserted, and the parameters vp and wa are omitted 
for the sake of simplicity. There is an obvious U(l) redun- 
dancy of this Hamiltonian. Let us denote the nodal Dirac 
quasi-particles field by Xq(x). The h/s gauge transfor- 
mation X xe'^T<P^ (f)^ + 2A, -> -f d^k leaves 
Eq. psp invariant provided that A is time independent. 
We then make the field redefinition from Xa to Xa as 
Xi — xie**^'/^ and X2 = X2e~"^ so that the nodal 
field becomes neutral under h/s gauge transformations, 
a 'nodon'j^ The above redefinition leads to a more con- 
venient form of the nodal Hamiltonian: 

" \ id+ -a+- ao) 

= -a+ + + ao)o-3 - i(?_(T2, (39) 

where the h/s gauge invariant field = — ^(9^^'* is 
introduced. 

Rotating the coordinate by 7r/2 successively, one may 
get the nodal Hamiltonian in the other three quadrants. 

D. Effective Action of 

In this subsection, we turn to the path-integral for- 
malism and derive an effective action (needed to discuss 
RVB gap equation) for in the nodal approximation 



by integrating out the holon fields. In the 1st quadrant, 
the effective Lagrangian in the Euclidean space for nodal 
quasi-particles is given by: 

=x{x)[rid^-2bf)]x{x), (40) 

where 7'' = {cr^:, -cr.y, ct^}, 9^ = {do,d+,d-} and bjf^ = 
{— ia+,Mo,0}. The effective action for 6^ (at T = 0) is 
defined as 

Si#K] = -lndet[7^(9^-z&f )]« 

J d'k J du;[bfultb'f]ik,u;). (41) 

By adapting the calculations of Ref. [s^, the leading terms 
of the bubbles for small lu, \k\,Lu/\k\ behave like 

Ull'^c,\k\, U^f+^c, Ulf^O. (42) 

The effective action in the other three quadrants is 
similar to that in Eqs. (|4ip and (|42|) . For example, the 
3rd quadrant can be obtained by rotating the coordinate 
by TT, therefore by changing a± — >■ —a± and d± — —d±, 
we can obtain the corresponding bubble H^i^ and gauge 
field bfj_. Note that the coordinate transformation does 
not involve the time axis. Then we have 

• H3'^d^Hi^;,if^ = 0,l,2, 

. bl''' = {ia+,iao,0}. 

Similar procedure can be applied to the 2nd and 4th 
quadrants. 

After summing over all four quadrants, in the 
quadratic approximation, we have the effective action 

«=+.- 

[a^Hooai + aoUuao - 2aoHoiai](fc,a;).(43) 

Using Eq. (|^^ and its analogues one sees that Eq. 
is a variant of the effective action for QED3. 

IV. SPINON RVB PAIRING 
A. Mean Field Lagrangian with Spinon Pairing 

In this subsection we derive the mean field Lagrangian 
for spinon pairing in the presence of holon pairs. 

In the PG region, the spinon part in the t ~ J model 
can be described by the massive sigma model in CP^ 
form Eq. (|17p. and the four fermion interaction term 
~ J{h\hyK'll){hjKKlj) (see the last term in Eq. pH)) ). 
is simply neglected for small doping 5 in considering the 
normal state, because it is proportional to 5^ . Note this 
interaction term is positive (for J > 0), hence repulsive 
due to the semionic mean field approach, contrary to 
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the usual fermionic case. However, once the holon pair- 
ing is stabilized, the gauge interaction between holon and 
spinon, overcoming the above repulsion forces the spinous 
to form singlet-RVB pairs and the above term becomes 
relevant. To investigate the spinon pairing, one can ap- 
ply a Hubbard-Stratonovich transformation to the four 
fcrmion interaction term, obtaining 

E + + h-c, (44) 

iv) 

where t = \{hihj) \ and in MFA 

At^ = ^r^{e'^^z,^z,p)^^r'm). (45) 

In the continuum limit we get the Lagrangian for spinon 
with a singlet spinon pairing 

/i=0,l,2 

I 



with = 0'* — 0" and = — ^d^cj)'^, both of which 
are h/s gauge invariant. The gradient of the </> field ac- 
tually describes the potential of standard magnetic vor- 
tices, since from Eqs. (|20l) and (|45|) (f) is the phase of the 
condensate of hole-pairs. 

By neglecting the gauge fields, one can work out the 
spinon spectrum, which can be obtained from the zeros 
of the determinant of the kernel T ^ in the momentum 
space: 

{^u^ + e + mlf-A Al,A^^%hk,^0. (48) 

We assume the rotational invariance for the spinon spec- 
trum, which requires 

Af^oA- +A-A5^o = 2<5,,|Af,op. (49) 

We can take Af g = Ag and A| q = ±iAg, where Ag 
can be a priori any complex number, and both plus and 
minus signs are allowed for A|g. Looking at the hole- 
pair order parameter Eq. (|20p . we see, however, that for 
consistency we have to choose the constant phases of Ag 
equal to —Bij to ensure the correct symmetry, being the 
hole-pair d-wave. From Eq. (|48p we obtain the spectrum 
for spinon: it has two (positive) branches 

E±{k) = ^JP +ml±2A'o\k\ (50) 
which arc plotted in Fig. 01 The positive branches of 



+ Y Ar{x)e"Pz^{x)d^zp{x) + h.c., (46) 

i=l,2 



where the index i in A''* labels the spatial directions 
and we set g and Vs to 1 for convenience. (The spa- 
tial derivative term in the square brackets has an im- 
plicit '-' sign, see Eq. P7)) .) As for the holon case, 
one can rewrite approximately the spinon pairing as 
A|(x) = Ai qc*'^''-^-' where is the phase of the spinon 
pairing amplitude. The Lagrangian Eq. (|46|) is invari- 
ant under the h/s gauge transformation Za — )■ ZqG*'^, 
Af^ — >■ Af^ + df^A and -I- 2A. It is not convenient 

to deal with the off-diagonal terms in the Lagrangian Cs , 
hence we transform the spinon field from Za to Za as 
zi = zie*'^°/^,52 = 22^"*'^°^^ so that the spinon field be- 
comes neutral under h/s gauge transformations. In terms 
of the new fields Za , the spinon Lagrangian can be writ- 
ten in a diagonal form ^s{x) = z''{x)Ts{x)z{x) where 
the 2x2 kernel reads (with Ag g = 0) 



(47) 

I 




I 1 1 ' 1 
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k 



FIG. 4: The spinon spectrum for S = 0.1. 



the dispersion are similar to those found in a plasma of 
rclativistic fermions;^ which suggests the following in- 
terpretation: if I A* I ^ the spinon system contains a 
gas of RVB spinon pairs, an analogue of Coulomb neu- 
tral pairs in the rclativistic plasma, either in the plasma 
phase if (A'') = 0, or in a condensate if (A'') ^ 0. For 
a finite density of spinon pairs there are two (positive 
energy) excitations, with different energies, but the same 
spin and momenta. They are given, e.g., by creating a 
spinon up and by destructing a spinon down in one of 



E -[^M-^(aA'+ 2^A<</')'^. -»MA^,g)tT, -iRe(A^_o)(T^]2-f m2-|A^ o|2, 

M=0,l,2 
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the RVB pairs. Notice that the minimum at J\k\ = |A*| 
in the lower branch is similar to the roton minimum in 
superfluid helium and has an energy lower than m^.; it 
implies a backflow of the gas of spinon-pairs dressing the 
"bare spinon" . Hence RVB condensation would lower 
the spinon kinetic energy. However, to make it occur 
one needs the gauge contribution to overcome the spinon 
repulsion generated by the Heisenberg term. 

B. Effective action of Gauge fields 

In this subsection, we derive the low-energy effective 
action of and (f> by integrating over the spinon fields. 
For this purpose we introduce a fictitious SU(2) gauge 
fields Yfj^ as follows: 

with 

/ ao + ao0\ 

y; = 2 Im(Alo) Re(Alo) a, + \d,^ . (52) 

Vlm(A|,o) Re(Ai,o) a2 + \d2<j^) 

Then the kernel Tg (see Eq. (|47p ) can be written in a 
compact form 

T,^Y.^^,^^Y;:^f+M\ (53) 

where we introduce the notation 

M = ^rn?, - 2|Ag|2 (54) 

for convenience. After integrating over the spinon fields 
Zoi, one obtains the effective action for and 9^0, 

5|ff[5^0,a^, A^] = lndet(r,) - ^1^, (55) 

where the constant term comes from the Hubbard- 
Stratonovich transformation. Since Eq. (|53p is formally 
describing a relativistic 2-component boson of mass M 
minimally coupled to the SU(2) gauge field F^, the lead- 
ing gauge invariant term is the Yang-Mills Lagrangian, 
i.e., the traced square (X]^i/ ^ly^/Ti/) of the field strength 
Y 

Y,. = y [a^r; - + s-'^y^y,'] (56) 

and one easily computes, with i = 1,2: 

yo---(ao + ao</')Re(A?o), 

^1^2 = Rc(A-Jo)(a2 +d2(f>)^ Rc(A^o)(ai 

Yy, = {ao + do(b)lm{Ato), 

Y,l = Im(A^o)(ai + 9x0) - Im(A-Jo)(a2 + 320), 



y^i = do{ai + di(j)) - di{ao + 5o0), 
2/i2 = 9i{a2 + d24>) - d2{ai + di(j)), 
+Im(A«io)Re(A^o) - Im(A^o)Re(A^o). 

Besides the Yang-Mills action there are also gauge non- 
invariant terms which arise from the ultra-violet diver- 
gences of the continuous model and must be included 
since the x, y components of Y^ are actually constant. 
For the 0th and 2nd order terms in and (f) we finally 
get 

s:^ - + Ein[(^' + Eimiu:^ + Elm, 

-t-|Agn2(ao + Idod^r + ia+ ^V^)^]}, (57) 

where a surface term (~ dia2 — 820,1) has been discarded. 
For I Ag| ^ 0, S'^V is the action of a gauged XY or Stueck- 
elberg model and the term in the last square bracket is 
the celebrated Anderson-Higgs mass term. 

C. Gap equation of spinon pairing 

The gap equation is determined by the saddle point of 
Ag] = 5:g°[Ag] + S:i[a, Ag] + ^,\[a] with respect 
to Aq. Note that since the interaction between spinous 
is repulsive, it is crucial to take the gauge fiuctuation 
S^g into account, unlike in the traditional BCS theory, 
where the electron interaction is attractive. To establish 
the gap equation for the modulus of the order parameter, 
we assume, as discussed e.g. in Refs. [33l - l35l for fermions, 
that one should neglect the phase (0) fluctuations. Let us 
also neglect for simplicity at first the holon contribution 
S^ff, then the resulting gauge partition function, denoted 
by Zg, is given by: 

^9 - ^K(-aV'^ + 9''9^-f |Ag|2Ana.](58) 

where g'^" = diag(l, 1, 1) and A'^'^ diag(l, 1/2, 1/2) are 
3x3 diagonal matrices and a cutoff A in both momenta 
and energy is understood. The equation of motion of 
gauge field reads 

-92a^ + 9'^(9.a.) + A'^''a, = 0, (59) 

which implies that (without source) satisfies the equa- 
tion 

A^-^a^a, = 0. (60) 

Eq. (|60p is a constraint, meaning that the massive vector 
bosons in two dimensions have two (physical) polariza- 
tion modes. Therefore, the calculation of the partition 
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function of the vector boson is not as trivial as integrat- 
ing over Op directly, in fact one must take care to count 
only the physical degrees of freedom. The details of the 
evaluation of Zg is given in Appendix |Bj Here we only 
give the final result 



Z„ = 



n 




(61) 



The first factor in r.h.s. of Eq. (pT|) contributes a constant 
term to the free energy which is neglected in the standard 
cases. However, in the present case it contains the spinon 
order parameter Ag, which can affect the total energy 
as varies, hence should be kept. Actually, SttM is 
the renormalization factor of the amplitude of the gauge 
action, which can be absorbed into by rescaling — 
V37rMap with a Jacobian left for the measure D\a^ in 



the energy-momentum space 



(62) 



is crucial, without which the gap equation has no solu- 
tion, since the last term is negative. In Eq. only the 
value of T is unknown, i.e., the nearest neighbor holon 
pairing strength, which is a very short range correlation 
and may not be accurate if being calculated via the long 
wavelength pairing A| in momentum space. However, we 

have already seen that extrapolating Aji to lattice scale, 
one gets the correct symmetry in real space. Hence we 
take Aq (fci^) as the value of r up to a scale factor. Let us 
briefly comment on the relation of our RVB gap equation 
Eq. ((64|) with that of the slave boson approach. Whereas 
in the slave-boson approach the RVB pairs are made of 
fcrmions and the Hcisenbcrg term is attractive, so the 
pair-formation is BCS-likc, in our approach the RVB 
pairs are made of bosons, and the Heisenberg term is 
repulsive, so the pair formation arises from the decrease 
in the free energy of spinous, via the lowering of their 
mass gap, induced by holon-pairing through the gauge 
field. 

So far we have not considered the vector boson quasi- 
particlcs, whose spectrum has two branches as derived 
from Eq. ([6T|) . 



where the power index 3/2 is due to the fact that is a 
three-dimensional vector. 

In fact in the spinon gap equation the term SttA/ al- 
ready balances the repulsive interaction. The contribu- 
tions of the spectrum of gauge quasiparticles, i.e.. the 
second and third terms in Eq. (|6T|) . do not change the 
spinon gap equation qualitatively. Therefore, for sim- 
plicity we focus only on the A'/-term, and the free energy 
including the contribution from the h/s gauge fiuctuation 
reads 



^^ln[(^2 + i?2(fc))(c.2+i52(fc))] 



3A3 

4 



, 2 2|Ag|2 
In m: ■ — 



Jr2 ■ 



(63) 



and contributes to the gap equation with the following 
term 



- y 

2V ^ 



1 



k,n=1.1 ^9 ^d,nn 2T 



(66) 



This contribution is positive and in balancing the gap 
equation Eq. ([M)) plays a role similar to the M-term. that 
turns out to be dominant. It is interesting to note that 
Eq. ((66)) is well defined in the gapped region Ag ^ 0, and 
if Ag = 0, it is proportional ~ TlnL which is divergent 
unless T = if the typical length of the sample L goes 
to infinity. Such an infared divergence seems to imply a 
first order phase transition when spinous begin to pair. 
However, this is not the case. In fact, when we take into 
account the contribution of holons S'^g to the action of 



gauge field (sec Eq. (j43 
(see Appendix [B]) : 



the dispersions (|64|) become 



It is straightforward to obtain the gap equation by taking 
derivative of Fg with respect to | Ag | ^ : 







3A3 
2^ 



A^ 
Jt2 



2|Ag|y 



E 



tanh 



2T 



i?+ tanh ^ J 



(.64) 



The first term originates from the gauge action due to 
the lowering of the spinon mass (m^ — ^ (m^ — |A''|2)i/2)^ 
while the second term comes from the original repulsive 
Heisenberg term and the last two terms are due to the 
spinon excitations. The first term in r.h.s. of Eq. (|64p 



i?W(fc) = WF + /(fc) + 



(fc) 




where C2 = 3ttMc2 and f{k) = SttMci^ Vprx,^ 

and the divergence disappears. 

In the low doping limit at T = 0, expanding the last 
terms in the r.h.s. of Eq. (|64|) we get 



A3/2 



1 



Jt2- 



(68) 
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As the doping 6 is decreased, r goes to zero faster 
than nis, because the spinon mass ~ |(51n(5| and 
t2 ^ ,5e~'=™''' (see Eq. which imphes that |Ag| has 

no nonzero solution for sufficiently small doping. In other 
words, there is a critical doping Sc at zero temperature, 
below which spinon pairing Ag must vanish. As the non- 
vanishing of Aq is a pre-condition for SC, this implies a 
critical doping for SC at T = 0. On the other hand, at 
the qualitative level, due to the cancellation of i5 between 
TOs and T^, if T ( i.e. the holon-pairs density) is suffi- 
ciently large Eq. (|68| does have a solution, because the 
remaining |lni5| is a decreasing function. Notice again 
the crucial role of this logarithm, coming from the long- 
range tail of spin-vortices. 

At finite temperatures, we need to solve Eq. (|64| nu- 
merically. The crossover temperature at which in mean 
field approximation Ag becomes non- vanishing is denoted 
by Tps (not yet the SC Tc) and is related to the for- 
mation of a finite density of RVB spinon pairs. From 
Eq. ([M)) we see that to have solution for the gap equa- 
tion we need r = (hihj) ~ Ag 7^ 0, consistently with 
the physical mechanism proposed, hence Tph > Tpg and 
when the spinon RVB pairs are formed together with the 
already formed holon pairs, producing a finite density of 
preformed hole pairs. Due to the (j) phase fiuctuations, 
however, although the modulus of the SC order param- 
eter A'^ ^ A" /A'' of (PD)) is non- vanishing, if the hole 
pairs are not condensed one cannot interpret it as the hole 
gap. The temperature dependence of As is presented in 
Fig. [Da. One can see that, although near Tps the behav- 
ior is the typical square root of mean-field, at low T it is 
definitely not BCS-like, never approaching a constant. 



V. SUPERCONDUCTIVITY 

Now we are ready to finally discuss the true SC tran- 
sition. 



A. Nernst crossover 

In this subsection we first consider the physical effects 
due to a finite density of hole pairs before their conden- 
sation. 

The gauged XY or Stueckelberg model of Eq. ([57)) is 
well known to have in the lattice two phases (see Ref. [s^ 
for a rigorous discussion, while Ref. for a numeri- 
cal analysis): Coulomb and Higgs. If the coefficient, 
^ |AqP of the Anderson-Higgs mass term for a is suf- 
ficiently small, the phase field cj) fluctuates so strongly 
that it does not produce a mass gap for and (e*"^) = 
in the Coulomb gauge (a gauge-fixing is necessary due 
to the Elitzur theorerr4i). This is the Coulomb phase, 
where a plasma of magnetic vortices- antivortices appears. 
In the presence of a temperature gradient a perpendicu- 
lar external magnetic field induces an unbalance between 
vortices and antivortices, giving rise to a Nernst signal, 
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FIG. 5: (Color online) (a) is the T — 5 phase diagram of 
the mean field gap equation of spinon for different values of 
MF spinon pairing A" (gray lines) which could be compared 
with different levels of the Nernst signal^i^; A" = is Tps. 
(The curves at high dopings are not quantitatively reliable as 
they do not take into account the crossover to the "strange 
metal"). The dashed line is Tph, the "upper PG crossover 
temperature". The dotted line is the crossover temperature 
between the pseudogap and strange metal phases, T* . (b) is 
the as a function of temperature for fixed dopings. The 
temperature and A" are in units of J. 



even if the hole-pairs are not condensed yet. Therefore 
we conjecture that this phase of the model corresponds to 
the region in the phase diagram of underdoped cuprates 
characterized by a non-SC Nernst signal and a compari- 
son between the experimental phase diagram in Refs. 00 
and the one derived in our model, supports this idea. The 
result is shown in Fig. [SJ where the thick lines are equal- 
Aq lines. One expects that the level of Ag is roughly 
proportional to the intensity of the Nernst signal and a 
comparison of the figure with the experimental data^ 
shows a qualitative agreement for the 5 — T dependence. 
Note that the Nernst data are strongly supported by the 



measured magnetic-field induced diamagnetic signal^ as 
well as by STM visualized pair formation^ and quasi- 
particle fingerprint4i The Tph line in the figure is the 
upper pseudo-gap crossover temperature determined by 
Aq(A:f) of Eq. (p4)) . hence it does not take into account 
the transition to the SM phase, therefore can only be 
taken as a qualitative trend. At extremely low doping 
{S < 0.03) the lines are not reliable because the quenched 
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approximation for vortices used in our approach is not 
valid for too fow vortex density. 
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FIG. 6: (Color online) The energy of the magnetic resonance 
Er estimated by 2ms J for different dopings, compared with 
the scaled critical temperature 3.2Tc. To compare with ex- 
periments we take J — lOOmeV. The inset is the experimen- 
tal results taken from Ref. 0, where the black solid line is a 
parabolic approximation to Tc rescaled by 5.8. 



B. The superconducting transition 

Now we consider the true SC transition. For a suffi- 
ciently large coefficient jA^' p, the gauged XY or Stucckel- 
berg model of Eq. (|57p is in the broken symmetry phase: 
the fluctuations of 4> are exponentially suppressed and 
(e"^) ^ at r = or there is a quasi-condensation 
(power-law-decaying order parameter) at T > 0; accord- 
ingly magnetic vortex-antivortex pairs become small and 
dilute, so the gauge field is gapped. At the same time the 
holon, and hence the hole, acquires the nodal gap, i.e. 
the gap outside the nodes. In fact, one can prove that, 
due to the fluctuations of the field cj)'^, in our approach a 
gapless gauge field is inconsistent with the coherence of 
holon pairs in PG, i.e., coherent holon pairs cannot coex- 
ist with incoherent spinon pairs, as sketched in Appendix 
C. On the other hand, due to the QED-like structure of 
holons-gauge action, the gauge field cannot be gapped (in 
all components) by condensation of holon pairs alone as 
shown by Eq. (|35]); only the condensation of RVB spinon 
pairs at the same time can open a gap to the gauge fluc- 
tuations and then the nodal hole gap. Thus as soon as 
e*"^ (quasi-)condenses, the same occurs to {hihj), so that 
SC emerges, since from Eqs. and ([57)) the SC order 
parameter is A'^ ^ A'^ / (hihj) ^ (Ag/Ao)e*'^ and now its 
modulus and the expectation of its phase are nonzero (at 



T = 0, or power-law decaying at T > 0). It follows that 

According to the above considerations, if we assume 
that the holon contribution to the gauge held is sub- 
dominant, as expected, the SC transition from the PG 
phase should occur roughly at a value of Aq determined 
by the gauged XY model. Then one can extract an es- 
timate of the critical value of Ag from a formula pre- 
sented in Ref. 14^ for the critical value of the coefficient 
of Anderson-Higgs mass term in Eq. (|57p . If one rescales 
the gauge field to have the standard coefficient 1/2 for 
the Maxwell term, denote by q the charge of the (j) field 
w.r.t. the rescaled and denote by /3 the coefficient of 
the Anderson-Higgs mass term, then such formula reads: 



(69) 



where /3c denotes the critical value. The value q = cor- 
responds to a pure XY model; in our case (see Eq. ([57]) ) 
q = 2v37rM. Unfortunately our Anderson-Higgs mass 
term is not isotropic in space-time, therefore to apply 
Eq. (j69p we symmetrize it, and a posteriori the precise 
choice of the coefficient turns out to be almost irrelevant; 
we choose (3 — | Agp/(127rA/). With this choice the so- 
lution of Eq. (|69p gives 



1287r2 



(70) 



and the choice of the symmetrized coefficient changes 
only the second almost irrelevant term. According to 
Fig. El one obtains for the SC state at T = a range of 
dopings from 6 w 0.04 to (5 ss 0.25. Tentatively extend- 
ing the formula Eq. ([70]) to finite T, one obtains for the 
critical temperature Tc the red dashed line in Fig. [5j For 
the critical value of Ag the value of M is quite small, 
hence q almost vanishes and within this approximation 
the SC transition is essentially of XY-type. This implies 
also that the gauge contributions of holons which have 
been neglected above would be actually self-consistently 
strongly suppressed. In general one can see from Eq. (|57p 
that in our approach a reduction of M, and hence of 
spinon kinetic energy, implies a reduction of the gauge 
fluctuations. The scale of in our approach is reduced 
w.r.t. a naive BCS value ~ Aj because a price has to be 
paid to overcome the spinon repulsion, so that its scale 
is essentially set by Aq. 

Let us now outline some physical consequences of our 
approach to SC that are presently under further investi- 
gation: 

1) In the SC state the gauge gap destroys the Reizer 
singularity (see Eq. ^T9\\ ) which is responsible for the 
anomalous T-dependent life-time of the magnon and elec- 
tron resonances in the PG normal state. Hence these res- 
onances become sharper at the SC transition. In turn, 
this improves the kinetic energy of the hole. There- 
fore in our approach the SC transition from PG is "ki- 
netic energy driven" as opposed to the standard BCS 
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"potential energy driven". The above feature is sup- 
ported by some experiments on optical conductivit y ^^'^^ 
and, in particular, a recent experiment on underdoped 
cuprateS)^ where one finds an increase of the kinetic en- 
ergy in PG, being consistent with our approach (due to 
the partial gap induced by the tt-Aux), and its sharp de- 
crease in the SC phase. This shows that within this gauge 
approach the compositeness of the hole, with a gauge glu- 
ing force coming from the single-occupancy constraint, 
proved to be essential in interpreting the transport^ and 
thermodynamicali^ properties of cuprate superconduc- 
tors, turns out to be a key feature also for the SC tran- 
sition. 

2) The appearance of two positive branches in 
the spinon dispersion relation for a suitable spinon- 
antispinon attraction mediated by gauge fluctuations (in 
particular those corresponding to the Z2 subgroup left 
unbroken by the condensation of the SC pairs) induces a 
similar structure for the magnon dispersion around the 
AF wave vector^ reminiscent of the hour-glass shape 
of spectrum found in neutron experiments 4^ Further- 
more, since the energy of the resonance is approximately 
twice the spinon gap Jtos ^ J(l — 25)|i51n5|^/^, it has 
a maximum in S near the maximum of Aq, and through 
Eqs. (inni) and ^U^, it is naturally related to T^. This 
appears as a key feature of our approach: because of 
the intimate relation between short-range AFO and the 
SC attraction, both coming from the same term in the 
representation of the t — J model Eq. ([TU)) (the third 
term, see also (|14p ). there is an intrinsic relation be- 
tween the energy of the magnon resonance and Tc- This 
feature qualitatively agrees with experiments;^ as shown 
in Fig. [51 



VI. DISCUSSIONS AND CONCLUSIONS 

Before concluding, let us briefly comment on the com- 
parison of the present proposal with other models on SC 
mechanism in cuprates. 

It is clear that our proposal differs in an essential way 
from the traditional BCS-Eliashberg approachf^ no mat- 
ter whether the electron-phonon interaction or the AF 
fluctuations serve as the pairing glue, SC being there "po- 
tential energy driven" . SC arises in our approach from 
PG exhibiting characteristic features of a doped Mott in- 
sulator, such as small FS, hence from the physical point 
of view this approach is an implementation of the basic 
ideas advocated by P.W. Anderson, attributing SC to the 
strong correlation effects in doped Mott insulatorsi^"— 
Furthermore, in our approach the leading part of the orig- 
inal Heisenberg term is used to provide the AF action 
for the spinous, by using the identity Eq. (holding 
for the bosonic spinons). Only the subleading term pro- 
portional to the holon-pair density is used to obtain the 
formation of a finite density of RVB-pairs in Eq. (|64p . so 
the derived SC can be viewed as vaguely reminiscent of 
Laughlin's gossamer SCi^ 



Our formalism shares some similarities with other ap- 
proaches, exploring the same underlying physical idea, 
with, however, some substantial differences. Both in the 
standard slave-boson^ and in the bosonic-RVB phase- 
g^j.jj^g57j58 approaches the Nernst effect and SC occur due 
to Bose-Einstein condensation (BEG) of bosonic holons. 
Since BEG persists for arbitrary small density in these 
approaches both Nernst effect and SC at T = occur as 
soon as the long-range AFO disappears. The same also 
happens in the standard "preformed pair" approaches^ 
due to the persistence of condensation of pairs in the ex- 
treme BEG limit. Instead, in our approach the repulsive 
interaction between spinons prevents the appearance of 
the Nernst effect below a critical doping, and the hole 
pairing occurs only when the holon pair density is suffi- 
ciently large to "force" the RVB spinon pairing via gauge 
coupling, while an even higher doping at T = is neces- 
sary to get SC. Similar "critical" dopings also appear in 
the phase-fiuctuation approach of Rcf. 3^, of which the 
main physical difference from ours is in that approach 
nodes appear in the Nernst phase, whereas in ours a fi- 
nite FS still persists and nodes appear only in the SC 
phase. Also, in the new version^ of the bosonic RVB 
phase-string model at T = a finite interval opens up 
between the long-range AFO and SC, due to the com- 
pact nature of the gauge fields; in this region, however, 
holons and spinons are "condensed" in contrast to our 
approach. 

Another peculiar feature of the approach presented 
here, distinctive from other approaches is the appear- 
ance of three distinct crossovers related to the PG phe- 
nomenology: in our notations Tph, Tpg and T* . The high- 
est one in T is Tph (the presence of t' there is relevant) 
where holons start to pair reducing the spectral weight of 
the holo2^ and producing, e.g., a deviation from linear in- 
planc resistivity. A lower one, Tps where incoherent hole 
pairs are formed, mainly affecting the magnetic proper- 
tics since a finite FS still persists, e.g., giving rise to a 
boundary of the diamagnetic/Nernst signal. Finally we 
have the crossover line T* crossing Tpg in the phase dia- 
gram; it is due to the peculiar phenomenon of the optimal 
TT-flux occurring only in bipartite lattices and it is not di- 
rectly related to SC. It corresponds to a change in the 
holon dispersion and is characterized by complete sup- 
pression of the spectral weight for holes in the antinodal 
region. Below T* the effect of short-range AF fluctua- 
tions become stronger and their interplay with thermal 
diffusion induced by gauge fluctuations gives rise to the 
metal-insulator crossover and the inflection point of in- 
plane resistivity. Such composite structure of crossovers 
seems also to emerge from recent experiments on optical 
conductivity!^ 

The relation found between Tc and the energy of the 
magnetic resonance might suggest that perhaps in some 
form at least part of the mechanism for SC presented here 
can apply also to SC materials different from cuprates, 
but with strong interplay between SC and AF. One pos- 
sible candidate is the recently discovered iron-arsenic 
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superconductors;^ which show similar phase diagrams as 
cuprates. However, the parent compounds in those sys- 
tems are not insulators, but rather semi-metals. On the 
other hand, the freshly found new iron-selenic systems^ 
do have insulating states as reference, and we can expect 
similar behavior to occur there. 

To conclude, in this paper the spin-charge gauge ap- 
proach is applied to derive superconducting properties 
from the t — J model with single-occupancy constraint 
describing the Cu-0 plane of undcrdoped cupratc super- 
conductors with the following distinct features: 

1) The same model and the same set of approximations 
are used to consider both normal and superconducting 
state properties without extra assumptions. The physical 
implications of the theory are derived explicitly, and in 
its totality are consistent with experimental observations. 

2) The interplay of antiferromagnetism and supercon- 
ductivity is taken as the underlying physical foundation 
and is implemented systematically for both normal and 
superconducting states. The same super-exchange term 
is giving rise to antiferromagnetism in the leading or- 
der, and producing superconducting pairing in the sub- 
leading order. As a consequence, a universal relation be- 
tween the superconducting transition temperature and 
the magnetic resonance mode energy is derived, in con- 
sistency with experiments. 

3) An unusual three-step scenario for the appearance 
of superconductivity is proposed: At the higher crossover 
temperature the charge carriers (holons) start to form 
pairs and they affect the charge transport properties (de- 
viation from the linear temperature dependence of re- 
sistivity); at the lower crossover temperature incoherent 
(local) hole pairs are formed, and the derived pairing am- 
plitude as a function of temperature/doping is consistent 
with the Ncrnst, diamagnctism and STM data; the true 
superconducting transition is derived as "almost" of the 
classical 3D XY-type, with a phase diagram in agreement 
with experiments. 
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For the sake of simplicity, we omit temporarily the sub- 
script a. One can introduce a unitary matrix A, 



A ^ ^ 



/I 10 \ 
-110 
11 

V 1 -i/ 



which transforms the matrix to 



A'^Hj:A 



V 



— Vpk 








-A'l \ 

k 





— /i + Vpk 


Ai 

k 










fi — Vpk 










yU 


+ Vpk y 



provided that the holon pairing parameter A^ is p-wave 
like, i.e.. A'' - = — A|. The spectrum of quasiparticles 
consists of two decoupled branches. 



•'±,k 



(«i.fc±Ai)2 + |A||2. 



The free energy at temperature T then reads 
^ |^ln(l-f-e"T^) + ln(l + eT^) 



F 



-T 



(A2) 



(A3) 



i=±,k 



According to Hellman-Feynman theorem, we have the 
gap equation for order parameter A|;, 



dF 
5A¥ 



y ^tanhf^") 
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(A4) 
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If we assume ^ > 0, the branch with energy -(as given 



in Eq. ([55])) is lower and responsible for the low energy 
physics of p-wave pairing. The corresponding quasipar- 
ticle field reads 



^k = ^i^k+k) 



(A5) 



Appendix A: Diagonalization of Mean Field 
Hamiltonian of holon pairing 

We introduce a four components spinor field, j: = 
(a r,b r,a^ -)*, in terms of which the holon 

V a,k' a,k' a-k' a.-k' ' 

Hamiltonian Hh a = 5^,-*^ -'H t with the 4 x 4 



In terms of "^-fields, the effective pairing Hamiltonian can 
be written as 

H':^ - Y^ivpk - ^^)^^J: - J(A^,V^1^1- + h.c), (A6) 



and the gap equation at temperature T can be obtained 
by neglecting the positive branch with i = + in Eq. (|A4p 
as already written in (P^ . 
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For the right Dirac cone, the p-wave pairing parameter 
takes the foUowing form in polar coordinate system 



Aj-== A(j(q) (cos 0,-- sine,-) 



(A7) 



with its angular and radial parts separated consistently 
with the gap equation. Substituting Eq. (|A7p into 
Eq. (|32p . we have at zero temperature 



A: 



8^ 



Jeff X 2kqcos{ 



(fc2 - g2)2 + ^ 2(A;2 ^ q2y^2 _^ ^^2^2 sij,2 



Ag (g) (cos Og' — COS 0,-) 



y^(/i - vpqy + [A'^ {q) (cos Off - COS 6^)] 



(A8) 



Note that the most important term comes from momen- 
tum around the Fermi surface, q ^ k and 6^ ^ 9j: or 
^5 ~ + ^fe J therefore we can neglect terms proportional 
to sin^(0j — 9q) in the denominator of the first fraction 
of the r.h.s. of Eq. (jASp . Then the angular part is just 
dropped off from the p-wave gap equation with only the 
radial part remaining 
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(AlO) 



where A is a momentum cutoff and E{x) and K{x) are 
the elliptic integral of first and second kink, respectively. 



Appendix B: Evaluation of tlie partition function of 
vector bosons 

In this appendix, we show how to compute the path 
integral of Eq. ([55]) of vector bosons including the contri- 
bution from the holon part (see Eq. p2)) ). The relevant 
Lagrangian can be decomposed into two parts 



,(-aV" + a^a,)a„ + a^m'^''a,] (Bl) 



where the holon's contribution is absorbed into mass 
term m'^'^ which in momentum space has the form 




with C2 = 37rMc2 and f{k) = ii^Mci^ v'j^k\ + v'^k^{see 
Eq. dSlD). 



It is not appropriate to integrate directly, because 
the redundancy of degree of freedom of vector bosons. To 
rule out the redundancy, we adapt a method developed 
by 't Hooft)^ We first rewrite the a-field in terms of the 
h/s gauge field A^^ and the phase field 0'': 



(B3) 



Clearly, is gauge invariant under the following h/s 
gauge transformation 



A 



An + dnA 



^ 0" + 2A 



(B4) 



The Lagrangian rewritten in terms of A^ and ((>''' is given 
by: 

SirMCg = A^i-d^g'"' + d''d'' +m'''')A^ 

-i^'^a^m^'-a.^" + 0'^m^'^a,,A, (B5) 

We choose the fixing gauge function as 

1 



(B6) 



whose derivative with respect to the infinitesimal gauge 
transformation Eq. (jBSp reads 



(B7) 



SF 

— = -m^^df^d, + l = D. 



Then, following the Fadeev-Popov-Dewitt's approach, 
the path integral involving only the physical degrees of 
freedom can be calculated as 



^9 = I m^A''] 



SF 



6A 



det[m°0(3^Af)3/2i?i/2^-i/2] 



(B8) 



with 



-a^g/^-^ + df'd" + m^'' - m^^'m-^^'a^'a^' (B9) 



K 



Note that in the determinant of Eq. (jBSp is orig- 
inated from the complete measure of £'[A*,Ae] of 
Hubbard-Stratonovich transformation. Finally, we ob- 
tain the following result 



(SttM)? 



t (w2 + fc2 + mii)^(w2 +^11 + m^k^)h 



-(BIO) 



which is given in Sec. IIVCI The poles of Eq. (|B10p lead 
to the spectra of gauge bosons in Eq. (|67)l . 



Appendix C: Coherence of holon pairs and gapless 
gauge field 

We sketch here the argument proving that the gapless 
transverse gauge field arising from Eq. (j42p is inconsistent 
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with the coherence of holon pairs in PG, i.e. with a 
non- vanishing expectation value of e^'^ at T = in the 
Coulomb gauge, implying that the (global) h/s symmetry 
is broken. Let us assume condensation of holon pairs, 
but not of spinon pairs, then the Anderson-Higgs "mass" 
term in Eq. ([57| at large scale simply renormalize the 
Maxwell term»^ Then in the Coulomb gauge the effective 
lagrangian for A^^ and cf)'^ has the following form: 



[co{A,iA + d^)A + AoAAo) 
-ciiAiV^A, + di(f>''VAdi(l)^) 

+C2(Ao + ao0'')'](^) (CI) 



+{c^' + {c„\k\^rY'^^W'hk,~co). (C2) 



Neglecting the subleading C2 term one can easily calcu- 
late the equal-time large-distance behaviour of the Green 
function of 6'^: 



G{x, = 0) = / d'^kduj 



Ci|fc|3 + Col/Cpw^ 



Therefore we have for large |a;| 



|x|i/^.(C3) 



with Ci,i = 0,1,2 suitable positive constants. Integrat- 
ing out the gauge field in the path-integral formalism we 
obtain the effective action for ip^ in momentum space 



''■](fc,c.)-0"(fc,c.)[ci|fc| = 



(C4) 



vanishing at large distance, so that the condensation can- 
not occur. 
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